Quantum circuit implementation of the Hamiltonian versions of Grover's algorithm 
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We analyze three different quantum search algorithms, the traditional Grover's algorithm, its 
continuous-time analogue by Hamiltonian evolution, and finally the quantum search by local adia- 
batic evolution. We show that they are closely related algorithms in the sense that they all perform 
a rotation, at a constant angular velocity, from a uniform superposition of all states to the solution 
state. This make it possible to implement the last two algorithms by Hamiltonian evolution on a 
conventional quantum circuit, while keeping the quadratic speedup of Grover's original algorithm. 
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I. INTRODUCTION 

While the standard paradigm of quantum computa- 
tion uses quantum gates (i. e., unitary operators) applied 
sequentially on a quantum register, recent developments 
have introduced a new type of quantum algorithms where 
the state of the quantum register evolves continuously in 
time under the action of some Hamiltonian. It includes, 
for instance, the "analog analogue" of Grover's algorithm 
or the quantum algorithms by adiabatic evolution that 
have been intensively studied lately 0i 13 ■ It has been 
shown that these Hamiltonian algorithms are genuinely 
quantum in the sense that they reproduce the quadratic 
speed-up of Grover's algorithm (see, in particular, the 
local adiabatic version of Grover's algorithm Hi). The 
purpose of this paper is, on one hand, to clarify the links 
between these Hamiltonian algorithms and their conven- 
tional discrete equivalents, and, on the other hand, to 
show how they can be implemented on a traditional quan- 
tum circuit. This second issue is important because it 
was never shown before how it could be done keeping 
the quadratic speed-up of Grover's algorithm. It appears 
that all these algorithms take a very similar form in the 
high dimension limit, which is particularly surprising for 
the case of the adiabatic search algorithm. Specifically, 
we see that the mixing parameter (which measures the 
mixing between the initial and final Hamiltonians in the 
adiabatic search algorithm) has to evolve in such a way 
that the instantaneous ground state rotates at a constant 
rate from the initial to the final ground state. This makes 
the link fully explicit with Grover's original algorithm. 



II. TRADITIONAL GROVER'S ALGORITHM 

First of all, let us briefly recall the principle of Grover's 
algorithm 0, • It is designed to solve the problem of 
finding the values x for which a function /(x), usually 
called the "oracle" , is equal to 1 (while it vanishes every- 
where else). As quantum gates have to be reversible, the 
quantum oracle must take the form: 



Of ■.nN®'H2-^n 



where the N candidate solutions 
states of the Hilbert space TLm, 



N ^'H2:\x) \y) \x) (gi \y ® f{x)) 

(1) 

\x) are taken as the basis 
while ® stands for the 



addition modulo 2. By considering the second register H2 
as an ancilla and preparing it in the state "^[|0) ~ 
the application of O / on both registers will be equivalent 
to the following unitary operation on the first one: 



Ur.nN^nN:\x)^{-iy^'''^\x) 



(2) 



To clarify the notations, we will throughout this article 
restrict ourselves to the case where there is only one solu- 
tion X = m (our results may easily be generalized to the 
case of M solutions, roughly speaking by replacing N by 
N/M in all the formulas below). In this case, /(m) = 1 
while f{x) = (Vx ^ m) and Uf may be rewritten 



Uf =I~2\m){m\. 



(3) 



Initially, we have no idea of what the solution could 
be, so we prepare the system in a uniform superposition 
of all possible solutions: 
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(4) 



This state may easily be obtained by applying an 
Hadamard transform H on the n = log2 N qubits re- 
alizing the quantum register Ti.N, initially prepared in 
state |0). The algorithm will also require the following 
operation: 



Uo = i/®"(/ - 2|0)(0|)H^" = / - 2|s)(s|. 



(5) 



Let us define as a Grover iteration the operator G — 
— UoUf. Throughout this article, we will always assume 
that 1 for simplicity reasons but also because the 

link between the different versions of the algorithm will 
appear much more clearly. In that limit, it may be shown 
that the Grover iteration becomes a simple rotation of 
angle in the subspace spanned by \s) and \m). More 
precisely, successive applications of G on the initial state 
|s) will progressively make this state rotate to the solu- 
tion state |m): 



G^ 



cos Is) -I- sin I m) (6) 
cosaf'^ls) -f sinaf'^lm) (7) 
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where 



^ dis 



2£ 

N 



(8) 



Thus, by applying the Grover iteration R'^^^ w jVN 
times, we obtain the sohition state |m) with a proba- 
bihty close to 1 with a quadratic speed-up with respect 
to a classical search, which would necessarily require a 
number of calls to the oracle f{x) of order N. 

Let us notice that as jVN is generally not an integer, 
we have to round it, for instance to the nearest lower 
integer, so that i?"^'^ = Lf V^J- This results in an error 



i'^l.)~|m)|l<sin 



that tends to zero for A'' ^ cxd. 
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(9) 



III. HAMILTONIAN EQUIVALENT OF THE 
ORACLE 



In the two Hamiltonian quantum search algorithms 
discussed below, we will use the Hamiltonian: 



Hf = I - \m){m\. 



(10) 



Let us show why it can be considered as equivalent to the 
oracle Uf. If we apply this Hamiltonian on a basis state 
I a:;) during a time t, it yields 



iHft 



x) 



\x) X = m 
e^**|x) Va; 7^ m. 



(11) 
(12) 



We immediately see that by taking t = n, we reproduce 
the operation —Uf, that is 



iHf 



-Uf 



(13) 



Conversely, it is possible to simulate the application of 
Hf during a time i by a quantum circuit using a one- 
qubit ancilla prepared in state |0), two calls to the oracle 
Of and an additional phase gate 

[/, = e-'*|0)(0| + |l)(l|. (14) 
Considering the circuit represented in Fig.^, we have 



Of[\x)(g>m - \x)(^\fix)) 



I(g>Ut[\x)(g>\fix))] 
Of [e-*(i--^("»*|x) (g> \f{x)) 



= e-*(i-/(=^))* 



x) ® \f(x)) 
® |0) 



= e-'"i'\x)®\Q). (15) 



\x) 
|0) 



— Of 



Ut 



Of — 



|0) 



FIG. 1: Circuit for implementing the evolution of a Hamil- 
tonian Hf during a time t by using twice the corresponding 
oracle Of. 



IV. ANALOG QUANTUM SEARCH 

Let us now consider the "analog" algorithm introduced 
by Farhi et al. in In addition to the oracle Hamilto- 
nian Hf, we will need a second Hamiltonian 

ffo =iJ«S"(/ - |0)(0|)iJ«®" =/ - (16) 

that is related to Uq as defined in Eq. in the same 
manner than Hf is related to Uf. The algorithm 
consists in preparing the system in the starting state 
l^an^^ = 0)) = \s) and then let it evolve under the time- 
independent Hamiltonian i7™ = Ho + Hf. We may show 
by simple calculation that: 



N 



(17) 

|m)] (18) 



e-'*[cosa""(i)|s) + jsina""(t)|m)], (19) 



where 



a^"(t) 



N 



(20) 



As in the traditional algorithm, the search works via a 
rotation from |s) to |m). However, the rotation is con- 
tinuous here, and follows a different path because of the 
presence of i in the second term of Eq. H19(l . The so- 
lution state is thus obtained with probability one if we 
apply i/'^" during a time T**" = ^ \/N. Let us also notice 
that 



a-(2j) 



(21) 



which shows that the application of iJ^" during a time 
rpan^j^dis _ 2 corresponds roughly to one Grover itera- 
tion. 

Suppose now we want to implement this analog algo- 
rithm on a quantum circuit. We showed in the previ- 
ous section how to reproduce the application of Hf with 
a circuit, but this does not allow us to directly apply 
JJ^^n = //^ _|_ ff j, Order to achieve this, we need 
to cut the evolution time T^'^ into i?''" small intervals 
AT — = -I such that we may approximate 



U{AT) 



by 



U', 



AT 



^-i{Ho+Hf)AT 



'iHoAT-iHfAT 



(22) 



(23) 
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Using the Campbell-Baker-HausdorfF approximation, 
which states that |||e^+^ - e^e^la G 0(|||[A, B]|||2) with 
|||A|||2 = max|||3.)j|^i denoting the operator norm 

of ^, we have 



|||C/(AT) - U'^A2 e 0{mo,Hf]\hAT^). (24) 



As 



the error introduced at each step is of order 



\IU{AT)-U'^rhe0{-^^ 



(25) 



(26) 



Since there are R^'^ successive steps, the global error 
made by this discretized analog algorithm is 



ll|C/(T)-([/;,^)«""|beO(^). (27) 



This simply results from the property that if the condi- 
tion jUj - U-iy < 1 is fulfilled Vj, then 

\ll[Uj-l[U',h<T.\lU,-U'^h- (28) 



We thus observe that for a number of steps i?''" ~ L"^J 
(of the same order in N as in Grover's traditional algo- 
rithm), we get the desired state with an error of order 
e. 

Furthermore, keeping the results of the previous sec- 
tion in mind, each step will have the same form 
g-iHoST^-iHfST ^ Grover iteration G = -UoUf, and 
therefore will require 2 calls to the oracle Of for being 
implemented with a quantum circuit. 



QUANTUM SEARCH BY LOCAL 
ADIABATIC EVOLUTION 



For this third algorithm, exposed in 3] and Q, we will 
once more need the Hamiltonians Hq and Hf, and we will 
initially prepare our system in a uniform superposition 
of all possible solutions \^'^^{t — 0)) = \s). This time, 
however, we apply a time-dependent Hamiltonian 



H{s) = (1 - s)Ho + sHf 



(29) 



to the system, where the mixing parameter s = s{t) is a 
monotonic function with s(0) = and s{T^) = 1. As |s) 
is the ground state of H{Q) — Hq, the Adiabatic Theorem 
tells us that during the evolution s{t), the system will 
stay near the instantaneous ground state of H{s) as long 
as the evolution of H{s) is slow enough. If this condition 
is satisfied, the system will thus end up in the ground 
state of H{1) = Hf, which is the solution state \m). 



Let us first study the path followed by \ip'^{t)) dur- 
ing the evolution. As H{s) only acts on the subspace 
spanned by \s) and |m) and as we start from \s), the 
path followed by |'0^'^(t)) will remain in this subspace so 
that the problem may be studied in this 2-dimcnsional 
Hilbert space. By calculating the eigenstates of H{s), we 
find 



where 



1^0 ;s) = 

\Ei;s) = 

Eo{s) = 

Ei{s) = 



N{Ei{s)-s)\s)+s\m) 
v/i?i(s)2 + (iV-l)(i?i(s)-s)2 

VNiEois) - s)\s) + s\m) 



N -1 

l-,/l-4^^s(l-,s) 



l-4^s(l-,s) 



(30) 
(31) 

(32) 
(33) 



If the adiabatic condition is realized (see 3), we have 

VN{Ei{s)~s)\s)+s\m) 



^E,{sr + {N-l){E,{s)-sy 



cosa''''(s)|s) + sina''''(s)|m) 



ad / 



where 



arctan 
1 



N{Eiis)~.s) 
2s 



■ arctan - , 

2 ^/N{1 - 2s) 



(34) 
(35) 

(36) 
(37) 



in the limit TV > 1. 

n/2 




FIG. 2: Rotation angle a^'^{s) for the adiabatic quantum 
search algorithm with A'^ — 32. 

The function q;^^(s) is plotted in Fig. [21 We see that 
the evolution is once again a rotation from \s) to |to), 
but which is not performed at a constant angular veloc- 
ity if s{t) is chosen to be linear in t, which corresponds 
to the quantum search by global adiabatic evolution orig- 
inally described in The observed angular velocity is 



4 



indeed greater for s close to 1/2 while it is smaller at 
the beginning and the end of the time evolution. Let 
us also notice that at discrete values s{t) = Sk, the con- 
tinuous path \ijj^'^{s)) coincides with the states of 
Grover's traditional algorithm. Thus, in the global adi- 
abatic search algorithm, the system exactly follows the 
path of Grover's algorithm, but at a varying rate. This 
suggests that this algorithm is not the correct adiabatic 
equivalent to Grover's algorithm. Moreover, we note that 
if s{t) — t/T^'^, then the Adiabatic Theorem imposes that 
rpad ^ 0{N), so that we loose the quadratic speed-up of 
Grover's algorithm (see 3). 

In order to circumvent this problem, we can perform 
a local adiabatic evolution as defined in Q . Then we get 
the solution state with an error less than e 

|||^ad(yad)^_|^^|| <^ (38) 

provided that we evolve at a rate such that 
N 



t{s) 



2e^/N - 1 
-t- arctan \/ N — 1] 



[arctan[VA^ - l(2s - 1)] 



N y/N^{2s - 1) + ViV^ 

arctan ■ 



1 - (iV- l)(2s- 1) 



2s 



arctan ■ , 

2e ViV(l - 2s) 



(39) 



in the limit N ^ I. Thus, for a local adiabatic evolution. 



et 

N 



(40) 



This is now a rotation at a constant rate during a time 
T*"^ — j^VN, so that we may consider this local adia- 
batic evolution as the right equivalent to Grover's algo- 
rithm. 

Let us study the implementation of this algorithm 
on a quantum circuit (we will closely follow the lines 
of the development exposed in Q). As for the ana- 
log quantum search, we discretize the evolution by cut- 
ting the time T^^ in R""^ intervals AT = Dur- 
ing each interval [tj^i,tj]{tj — jAT), we approximate 
the varying Hamiltonian H{s{t)) by the constant one 
Hj = (1 — Sj)Ho + SjHf {sj — s{tj)). In order to 
evaluate the error introduced by this approximation, we 
will use the following lemma: 

Lemma 1 Let H{t) and H'{t) be two time- dependent 
Hamiltonians for < t < T , and let U{T) and U'{T) 
he the respective unitary evolutions that they induce. If 
the difference between the Hamiltonians is limited by 
\H{t) — H' {t)\2 ^ <5(0; then the distance between the in- 
duced transformations is bounded by \U{T) — [/'(T)|||2 < 



^2]^mdt. 



This lemma is a straightforward generalization of the 
one introduced in 3], with the important difference that 



the Hamiltonian difference 5{t) may vary in time. This 
will be crucial in order to keep the quadratic speedup of 
Grover's search after this discretization procedure (this 
is reminiscent to the distinction between the global and 
the local adiabatic evolution) . The proof of Lemma 1 is 
left to the reader. 

Thus, the approximation we made is equivalent to 
replacing the actual Hamiltonian H{t) = H{s(t)) by 
H'{t) = H{s'{t)) where s'{t) is a new monotonic func- 
tion approaching s{t) but varying at times tj only (see 
Fig.EJ. 




FIG. 3: s(t) and its discrete approximation s'{t) (using R^"^ 
20 steps) for the local adiabatic algorithm with A'^ = 32. 



We have 
\\H{t)-H'{ 



h = l\H{sit))-His'mh (41) 
= W{t)-s{t)\ \mn~Hfl\2 (42) 

< {-AT-^0{AT'))\m,,^Hjh. 



We may use Lemma 1 with S{t) = ^ATfH^n - iJ/ib 



S{t)dt = AT\IH,„-Hf\h ^ 



dt (43) 



- AT\IH,,,,-H 



/II 2 



ds 



< AT 

as we may easily see that 
lemma gives: 



(44) 

(45) 
(46) 



Hfiy < 1. Now, the 



l\UiT) ~ U'iT)\h < V2AT = y 2-^ (47) 

so that in order to keep an error of constant order e for 
growing N, we must choose a number of steps propor- 
tionnal to T^^, that is R""^ = [^J. 

For each step, we have to apply Hj during a time AT = 
, that is the unitary operation: 



TTi _ -iHjAT _ -i(l-Sj)HoAT-isjHfAT 

— c — e 



(48) 
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Grover 


Analog 


Adiabatic 


#steps R 
Step j 

State \ipj} 
Angle Qj 

Error \\\->Pr) - \in)\\ 


Stoj — Stfj = IT 

COS aj\s} + sin aj\m) 

2j 


cos aj\s) + i sin aj j m) 

Hj TV 

ViV 2 
0(6) 


COS Qj |s) + sin oij\m) 

2 

0{e) 



TABLE I: Summary of the properties of the three quantum search algorithms. Here SIq^j and 5t fj are the times during which 
Ho and Hf have to be applied in step j. Notice that we have omitted some irrelevant global phases in front of and |V'-r) 



As for the analog algorithm, we use the Campbell-Baker- 
HausdorfF approximation and replace J7,' by 



(49) 



The error introduced at each step by this approximation 
will be 



lie/; ~c/;ill2e 0(5,(1-.,)-^). (50) 

For the R^'^ steps we have U'{T) = t^jj so that using 
Eq. |j2HI) gives 



(51) 



\\U'{T)-\{u'^h^oCj^, 



that have been found so far are very closely related. They 
all perform a rotation from the uniform superposition 
of all states to the solution state at a constant angu- 
lar velocity, even though a slightly different path is fol- 
lowed by the analog quantum search algorithm. Their 
similarities become even more obvious when they are 
implemented on a quantum circuit as they all require 
a number of steps of order ViV, each step having the 
same form e-iHoSta^-iH^stj ^ ^^^^ ^^isX the "duty cycle" 
Stf/{Sto + Stf) varies along the evolution according to 
a specific law in the case of the local adiabatic search 
algorithm. Finally, we have shown how one can realize 
these basic steps on a quantum circuit by using two calls 
to the quantum oracle. These results are summarized in 
Table H 



Consequently, the number of steps i?^"^ — [-^J required 
in the previous approximation (i.e., replacing H{t) by 
H'{t)) results in an error of order here. 
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